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Abstract 

Four stroke Otto refrigerator cycles with no classical analogue are studied. Extremely short cycle 
times with respect to the internal time scale of the working medium characterize these refrigerators. 
Therefore, these cycles are termed sudden. The sudden cycles are characterized by the stable limit 
cycle which is the invariant of the global cycle propagator. During their operation the states of the 
working medium possess significant coherence which is not erased in the equilibration segments due 
to the very short time allocated. This characteristic is reflected in a difference between the energy 
entropy and the Von Neumann entropy of the working medium. A classification scheme for sudden 
refrigerators is developed allowing simple approximations for the cooling power and coefficient of 
performance. 

PACS numbers: 03.65.Yz,05.70.Ln, 07.20.Pe,05.30.-d 



1 



I. INTRODUCTION 



Magnetic cooling was initiated more than a half century ago l). The cooling agent 
paramagnetic salts was cooled both by adiabatic demagnetization [l| and by adiabatic mag- 
netization Q depending on the sign of the isothermal change of entropy as a function of 
the magnetic field. Currently, adiabatic demagnetization is an efficient technique for cooling 
detectors in space missions and also for home refrigerators without freon gas 3-6]. These 
devices are realizations of quantum refrigerators a subject which has been of growing interest 
in the last decade jz-24]. In this study we explore an extreme mode of operation of mag- 
netic refrigerators where the short time allocation on the cycle segments leads to quantum 
characteristics. 

A prototype of a magnetic four stroke refrigerator is investigated. The working medium of 
;he refrigerator is a magnetic material modeled as an ensemble of two coupled spin systems 
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26| . The working cycle is composed of four segments: two isomagnets, one magnetization 



and one demagnetization segment. The dynamics is generated by a completely positive map 



271 ] . which settles in a self repeating cycle - the limit cycle [281 ]. 

In the microscopic model two timescales emerge, the cycle time and the internal timescale 
determined by frequencies of the working medium. When the allocated time on each segment 
of the cycle is large compared to the internal timescale, the cycle is quasi-adiabatic. Under 
these circumstances the state of the working medium is close to an equilibrium or Gibbs 
state. These states are characterized by the expectation value of the energy. As a result the 
state p is diagonal in the energy representation [p, H] = 0. 

When the cycle period becomes comparable to the internal timescale, the state of the 
system lag behind the change of the external parameters. The cause of this effect is that 
the external control Hamiltonian does not commute with the internal Hamiltonian. In this 
case, the states of the working fluid can not follow the instantaneous energy levels, therefore 
additional energy is stored in the working medium. This additional energy is accompanied 
by large off-diagonal elements of the state of the working medium p. The dissipation of this 
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additional energy is the quantum analogue of friction 

A quantitative distinction can be made between two sets of cycles depending if the time 
spent on each segment is shorter or longer than the internal timescale This distinction is 
a boundary range between the two sets. The sudden cycles are characterized by segment 
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times which are much shorter than the internal dynamical timescale. The regular cycles are 
those with segment times which are longer than the internal timescale. Our previous studies 



focused on regular cycles 
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29| . Intermediate cycles exists where some segments have a 



short time allocation compared to the internal timescale and other segments have a long 
time allocation. 

n the regular cycles the different segments can be characterized and studied individually. 
In 2^1 Fig. 4 is a representative example of those cycles where the different segments have 
many cycles, and can be treated almost independently, as opposed to the sudden cycles of the 
recent study. Typically, the states of the regular cycles are almost diagonal in the energy 
representation. This energy dominance is even more pronounced at the contact points 
between the cycle segments. The optimal cycles of the regular type are strictly diagonal 
in the energy representation on the contact points. These cycles are termed frictionless 
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The sudden cycles are characterized by large off-diagonal elements in the energy repre- 
sentation [p, H] 7^ 0. As a result, the sudden cycles are characterized by the Von Neumann 
entropy Svn = —tr{plnp} being different from the energy entropy Se = ~YlPn^ n Pn 
where Pi the expectation value of the nth energy level. Se > Svn with equality only when 
[p, H] = 0. These cycles are different from the time optimal bang bang type cycles which 
have significant off-diagonal elements in the energy representation on the adiabatic segments 



but are diagonal on the contact points 
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The large off diagonal elements of the density operator p in the sudden cycles accompany 
their global behavior. The elements do not vanish on the connecting points between seg- 
ments. As a result additional global links are generated. The interrelation of the observable 
values generate off diagonal elements of p. Increasing further the cycle time, the amplitude 
of the off-diagonal parts is suppressed. The coherence, defined in Sec. II characterizes the 
suddenness, showing that it decreases with increasing cycle time. 

Another important property of the sudden cycles is that their performance is not neces- 
sarily optimal. Starting from an arbitrary initial state the relaxation toward a limit cycle of 
sudden cycle might take several thousands of iterations, as opposed to the regular cycles, 
where generally several iterations are sufficient. When the sudden cycle is closed, and leads 
to refrigeration, there is a close neighborhood of parameters which lead to similar cycles. 
These neighborhoods of analytic behavior are very small. As a result there are islands of 
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parameters which lead to sudden refrigeration cycles which are disconnected from other 
islands. 

To gain insight on the sudden cycles, analytical approximate expressions for the cooling 
power of the cycle, its coefficient of performance (COP), the coherence as a function of 
time, and the entropy generation of a cycle are obtained It is shown, both numerically 
and by analytic approximations, that the cooling power of the sudden cycles achieve a 
maximum value as a function of the inverse temperature multiplied by the coefficient of 
the inner coupling, as opposed to regular cycles which depend exponentially on the inverse 
temperature. 

II. THE CYCLE OF OPERATION, THE QUANTUM HEAT PUMP 

A heat pump extracts heat from a cold reservoir, and transfers it to the hot reservoir. The 
operation of the heat pump is determined by the properties of the working medium and the 
coupling to the cold and hot baths. The cycle of operation is defined by the external controls 
which include the variation in time of the field with the periodic property u(t) = u>(t + r) 
where r is the total cycle time synchronized with the contact times on the different segments 
of the cycle. The cycle studied is composed of the following four segments, see Fig. |2j 

1. Segment A — )■ B (termed isomagnetic or isochore ), the field is maintained constant 
bj = uj c , the working medium is in contact with the cold bath of temperature T c 
with heat conductance T c , and dephasing parameter 7 C for a period of r c , and with 
propagator U c . 

2. Segment B — >■ C (termed magnetization or compression adiabat ), the field changes 
from u c to uih in a time period of r ch with propagator U ch . 

3. Segment C — )■ D isomagnetic , or isochore , the field is maintained constant u> = u>h, 
the working medium is in contact with the hot bath of temperature Th with heat 
conductance T h , and dephasing parameter ^ h for a period of and with propagator 
U h . 

4. Segment D — » A demagnetization or expansion adiabat , the field changes from u>h to 
uj c in a time period of r hc , with propagator Uh c - 



4 



In the basic paper, [101 ]. section III describes in detail the definitions of the segments 
isochore and adiabat . A somewhat different approach can be found in the paper of Quan 
et al, [18] . 



At the limit cycle, all the values of the expectation values exactly repeat themselves 
during each cycle time r. The propagator of the cycle will be termed U g i b a i, where U g i bai is 
constructed by the individual propagators (for example at point A in Fig. [2]), as: 



global 



(1) 



The limit cycle is characterized by an invariant eigenvector of Ugiobai , with eigenvalue lfone' . 



The dynamics of the refrigerator's working medium follow our previous studies [10|, [28|, |29| . 
We construct the segment propagators U by solving the dynamics for quantum thermody- 
namical observables. This dynamics is g enerated by completely positive maps within the 



formulation of quantum open systems 
£ in the Heisenberg picture, 



31 



321 ] . It is generated by the Liouville superoperator, 



dA 
~dt 



z[H,A]+£ D (A) + 



OA 
~dt 



(2) 



where Co is a generator of a completely positive dissipative Liouville superoperator, which 
includes the temperatures T c /h of the reservoirs. The state of the system p is then recon- 
structed from a finite set of observables. 

The Hamiltonian of the working fluid has the structure: 



H = Hint + H e xt(t) 



(3) 



where [Hj nt , ~H. ext {t)\ ^ 0. We chose a cycle that when the working fluid is in thermal contact 
with the reservoirs the Hamiltonian is stationary. 



Motivated by a refrigerator operating with magnetic salt 



3j-[6j, a simple model was con- 



structed with a working fluid composed of pairs of coupled spins. The uncontrolled, internal 



Hamiltonian becomes (h =1): H 



int 



) 



JB 2 



where <x 



represents the spin-Pauli operators, and J scales the strength of the inter-particle interac- 
tion, which is assumed to be constant, for a given pump. When J = 0, the model represents 
a working medium with noninteracting atoms jsj. The external Hamiltonian is chosen to 
be: H exi = \oj{t) (&l <8> I 2 + I 1 <g) <J 2 ^j = uit)^ where u(t) represents the external 
control field. The total Hamiltonian becomes: H = u(t)'B 1 + JB 2 



The eigenvalues of H are e% = -Q(t), e 2/3 = 0, e 4 = n(t) where Sl(t) = ^uj{t) 2 + J 2 , 
which is the temporary energy scale, which is dependent on both the internal and exter- 
nal parts of the Hamiltonian. At various times H(£) does not commute with itself since: 
[B 1 ,B 2 ]=2iB 3 ^0, B 3 = l(*l®a 2 x + al®&l) 



The explicit equation for the dissipative part Ld(A) [3l| is 

C D (A) = ~ ^(FjFjA + AF j F})') , (4) 

j 

where the (Fj) were identified with the raising and lowering operators, from energy level 
n to n — 1 and vice versa. These operators were constructed by first diagonalizing the 
Hamiltonian, then defining the (Fj) operators in the energy representation. The lowering 
transition rates n , h were chosen to be equal for all the four transitions, while the raising 
transition rates K , h were obtained by forcing detailed balance relation on the hot and cold 
isomagnetic segments i\./ h / i^, h = exp(—Q r /h/T r /h). The heat transfer rate to the baths 
will be T c /h = K c i h + ^ c j h i (&b = 1) lflf, see also the Appendix. 

In addition a dissipative generator of elastic encounters is added described as 

C D .(A) = - 7 [H,[H,A]] (5) 

where 7 is the dephasing constant. 

The state of the system p can be expanded by a complete set of orthogonal operators 
on the Hilbert space of the system. The expansion coefficients are proportional to the 
expectation values of these operators. This set can form a complete vector space to represent 
the propagators U{ on different segments. 

A thermodynamically inspired set of observables is a minimum set of operators which 
completely defines the state of the working medium when it reaches the limit cycle. For 
thermal equilibrium the energy and identity operators are sufficient. For the limit cycle this 
set has to be expanded. The set is initiated from the energy H and new operators are added 
which are dynamically coupled to the energy. This set is formed from a linear combinations 
of the stationary closed set |b| of operators: 

H = w(t)Bi + JB 2 L = - JBi + u(t)B 2 C = fi(t)B 3 (6) 

The three operators defined in Eq. form a closed Lie algebra, for they are linear combi- 
nations of the original operators Bi, which also form a closed Lie algebra. 



To uniquely define the diagonal part of the state p in the energy representation, the 
original set has to be supplemented with two additional operators: V = fiB 4 = ® 



& 2 Z - I <g) &l) and D 
expanded as: 



With this operator base the state p can be 



p = h + I ((H)H + (L)L + (C)C + (V)V + (D)D) (7) 

V and D commute with H. The equilibrium value of (V) is zero, and once it reaches 
equilibrium it does not change during the cycle dynamics. This means that on the limit 
cycle the state p can be reconstructed by four expectation values: E = (H), L = (L), 
C = (C) and D = (D). In the energy representation the state p becomes: 

fl + UD-2E) U L + tC ) \ 








o i-^o 

l-i£> 

\ ^(L-iC) 1 + ±(D + 2E)J 

A measure of the off diagonal elements in the energy frame is the coherence 



C = tr{(p e -p ed ) 2 } 



(8) 
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(9) 



where p ed is the diagonal stationary part of the density operator in the energy frame. ^From 



Eq. 



the coherence becomes C 



L 2 +C 2 
2C 2 



. Fig. |3] shows the transition from a sudden 



cycle to a regular one on the adiabatic segment. An increase in the time allocation on this 
segment decrease the coherence C. 

The vector space defining the propagators is constructed from the four operators 
H, L, C, D and the identity I I = (H, L, C,D,I). Using this set the propagator Ui(r) 
on the isochores (or equivalently on the isomagnets) becomes [101 ] : 
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-^f (e- r -/" T - 1) 



e (- r c/h T ) 






where K = e l -~ ]LTc / h+lc ' hn2 ^ T \ T c / h , j c / h , Q, r are defined above. 

The periodic functions in Eq. f llUp mean that the isomagnetic segments are quantized. 

27r, the two coupled equations of L, C complete a period. The other two 
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FIG. 1: (Color online) Top (a): An example of an extreme sudden cycle with analytical scheduling, 
in the (Q,Se) plane (curved shape ) and (Q,Svn) plane (shrinked rectangle shape), together 
with the isotherms corresponding to the cold/hot baths temperatures, T c /T^. The cooling power 
Q c /t = 1.2 • 10 -6 . The cyle parameters are: T c = 14, = 15, J = 2.,uj c = 0.1, uj^ = 6.. The time 
allocations: T c ,r ch ,T h ,T hc = 0.9,0.00035,0.00025,0.00035,4 = 0.36, 4 = 0.328, Bottom(/3): the 
corresponding cycle trajectory in the H, L, C space. 

expectation values E and D are decoupled from L, C the expectation values of the operators 
L and C. Quantization exists also on the adiabats . Closed form solutions for the adiabats 



leading to quantized motion are obtaind for a constant adiabatic parameter \i = j£ [29 1. 



Sudden cycles can be classified according to the time spent on the different branches. Fig. 
[T] displays the most extreme case. It seems that the two magnetization and demagnetization 
adiabats almost coincide in the entropy frequency plane. Also one recognizes that the cold 
isotherm almost touches the only point where the Se and Svn meet, showing that further 
cooling is impossible. The Von Neumann entropy Svn — ~ tr{p\np} is constant on the 
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FIG. 2: (Color online)A typical sudden cycle of the refrigerator with additional time allocation on 
the segments. Top(a): in the (Q, Se), and (O, Svn) plane (rectangle), together with the isotherms 
corresponding to the cold/hot baths temperatures, T c /Th. Note Se > Svn- The cycle parameters 
are: T c = 2.175, Th = 2.9, J = 2.5, uj c = 2.5, ujh = 10.. The time allocations are: T c ,T c h,Th-,Thc = 
0.2, 0.21, 0.44, 0.21, k^ c = 0.36, 0.328. Bottom(/3): the corresponding cycle trajectory in the H, L, C 
space. The longer time allocation on the hot isomagnetic segment and the larger frequency 
allows to complete approximately 3/4 of a period. 

adiabats since the dynamics on these segments is unitary. The bottom of Fig. [TJ shows the 
trajectories in the H, L, C space indicating that the dependence of the different segments 
forces the trajectory to reside on a plane. 

Fig. shows a second type of sudden cycle with addional time allocation on the isomag- 
netic segments. The bottom of Fig. [2] shows the trajectory of the cycle in the H, L, C space. 
One notices the global property of the cycle, by realizing that the isomagnetic segments 
complement each other (see also the third sudden cycle, Fig. |4]). The rotation of L and C 
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FIG. 3: (Color online) The coherence C, Eq. Q as a function of £1 for cycles with different time 
allocations on the adiabats . The cycle parameters correspond to the data of Fig. [T] except the 
different The = r c h- Notice the decrease in coherence when the time allocation on the adiabats 
increases. 

complete one period split between the hot and cold isomagnetic segments. In addition, Fig. 
[2] shows that further cooling is possible, because the end point B of the cold isomagnetic 
segment is below the point E where the vertical BE intersects with the cold isotherm. Point 
E is the upper-bound on Se for a cooling cycle. The maximum heat per cycle that still can 
be extracted from the cold reservoir is T c (Se(E) — Se(B)). 

One should also notice that in all the sudden cycles, the significant difference between 
gmax _ gmin anc j gmax _ ^mm^ or generally the difference between Se and its corresponding 
lower bound Svn-, all of which indicate the large off diagonal elements of the density matrix 
in the energy representation. The reason is that while H, L, C are created on the adiabats, 
the very short time on the hot isochores is not sufficient to equilibrate the energy of the 
working fluid. Generally, the time on the segments is so short, that less than one period is 
achieved. 

III. ANALYTICAL APPROXIMATIONS FOR THE COOLING POWER. 

The cooling power is the object of this study. The amount of heat extracted from the 
cold bath is defined by the difference in energy of the working medium between point B and 
point A Cf. Fig. The cooling power, V c = Q c / T i is the heat flow from the cold reservoir 
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FIG. 4: (Color online) A sudden cycle with additional time allocated on the isomagnetic segments. 
The cycle parameters are: J=2.5, = 2.9, T c = 2.175, = 10, oj c = 2.5 1= 3.35, k c J,= 0.328 
T h = 0.442, r hc = 0.00744, r c = 0.527, r ch = 0.00824. Top(a): In the (Sl,S E ), (^,S V n) planes. 
Bottom(jS): The corresponding trajectory in the H, L, C space. Linear scheduling of u(t). 

into the system divided by r. 

V c = QJr = (E B -E A )/r. (11) 

The values of Ea and Eb are calculated from the limit cycle invariant vector X of the cycle 
propagator U cyc at these points. 

The commutator of the cycles' propagators: 

[UabMc] (12) 

where Uab = UhMhUch supplies an indication of the cooling power. When \UabMc\ — 0, 
<2 C = 0. This relation is used to check approximations of the propagators U. The sudden 
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limit on the different segments suggest a short time approximation. We use the commutator 
Eq. f[T2T) to check the minimum order of the approximation. 



A. Approximations for the adiabats 



The exact solution for the propagator on the adiabats , U a d for a constant adiabatic 
parameter /x has been derived in Ref. |29|, Eq. (18): 
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(13) 



where g = a/ 1 + /i 2 , s = sin(qQ) and c = cos(g0). The angle for the hot-to-cold adiabat 



is defined as: 



©he = T adi (l/K hc ) (arcsin(o; c /fi c 



axcsm(u) h /Vt h )) = ^-(^ c 



A' 



fee 



(14) 



A similar expression is obtained for cold-to-hot angle Q c h, where: K^ c = j(u c /Q c —0Jh/ft<h) — 
- ~j(uh/^h ~ uj c /Q c ). The propagator Eq. ( 1T3|) is proportional to the compression 



ratio: C — jf- except for the last term which corresponds to the identity. 



Expanding the terms composing Uh c for short time (r c h = Th c 
= q. This result follows from the definitions of /i and q\ \x hc 



T a di < 1), one obtains 



Khc 



K ch 

1~adi 



where the relation between /i and r adi is obtained from [29j. Since r a di <C 1, from the 
definition of fi one obtains that ^> 1. Further: q = \Jl + /i 2 . Therefore = q. Under 
these conditions Uh c , Eq. ([TBI , is approximated as: 



(appr) 



he 



/ c -s \ 

s c 

10 

1 



(15) 



\o o 00|/ 

The explicit expression for the argument qQ becomes: qQ = |/i|0 = |i^|/r di = 

l-^l/'T'adiwfC^/ic) — (^/ic)- A similar expression is obtained for . Notice that in Eq. 
f[T5"|) . only the vector components H and L exchange their values. 
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1 . Further classification of sudden approximations on the adiabats 



A further classification of the propagators on the adiabats is defined according to the 
values of the operating parameters u c , Uh and J is based on additional simplification of 
either Eq. (TLj) or Eq. (fI5]b 



1. When co c <^ J, Uh ^> J and also r adi <C -Mx in the range of t, then the propagator 



(appr) 
he 



will simplify to the form: 



U, 
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he 
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(16) 



2. When u c ~ J, Uh ^> </ and also r ad j <C j^y. In this case will have the form : 

< 75 A ^ 
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3. case(a) 



When u c , Uh^> J and also T a di <C ^y. then ttj^ becomes: 

/ 1 \ 
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z^ a) » ^ 10 
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Case(b) 

When w c , J and also w 1. 

^From the condition w c , Uh^> J the ratio ^ becomes: 
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:i/2)(-) 2 . 



(17) 
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^From the definition of \i in Eq (1TB]) . K adi after Eq. (fT4"]) and condition (b), follows 
^adi < 1) therefore \i < 1. In addition from the definition of g follows g = 1. Therefore 
Eq. f|T3|) simplifies by neglecting /i 2 to the following propagator: 



(36) 
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(20) 



In ref. 



29], when expanding H,L,C in Eq. (TT5|) to first order in /i, a propagator 



similar to Eq. (1211 was obtained. 



^From g = 1, Eqs. f[T9"|) and (fT4l) . the argument gO of the trigonometric functions 
becomes 

(JV2) (l/aj - 



aai (J/2) - l/ul) ™" 

In this case the argument gO = of the trigonometric functions is large therefore 
it cannot be approximated to first (or second) order. 



B. Approximations for the isomagnetic segments . 



For the isomagnetic segments the propagators U c and Uh are approximated as follows. 



1. First we assume that the time allocation on both isomagnetic segments is short enough, 
so that in Eq. ( ITU]) , sin(fir) = 0, and cos(fir) = 1. In addition, 7 = 0. Then Eq. 
ffTO"]) simplifies to: 

/ e (-Tc/ h T) E eq {l - e { ~ r c/ hT) ) \ 

e (-r c/h r) 

U) -00 e (" r ^ r ) (|22) 

±(E eq (e- r °/>> T - e- 2r ^ r )) 



e -2r c/h r ^£^ e -r c/h r _ ^ 
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/ 



2. On the isomagnetic segments we assume short time on the hot isomagnetic segment, 
leading to a first order approximation in time. The result is an approximate propagator 
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for the hot isomagneticic segment: 

( l-T h r h 
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(23) 



1 

A first order approximation also on the cold segment leads to the commutation of the 
segment propagators Eq. (fT2~j) to vanish \UabMc\ — 0. Therefore for the cold segment, 
the unapproximated Eq. ffTOj) is employed. 

Another possibility is that the hot and cold isomagnetic segments are swapped, and 
the time on the cold isomagnetic segment is short, while the hot isomagnetic is not 
approximated. 

3. In this case both of the isomagnetic segments are approximated to the second order 
in time leading to: 

u c/h ~ 



( + f 
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1-2G + 2G 2 ^(G-f 



and G = Ft. 



C. The approximate cooling power 

The heat removed from the cold bath at each period is calculated from the global propa- 
gator, Eq. ([1]). The approximations of the segment propagators are employed, Eqs. fTT6l) to 
(|24|) . to obtain a global propagator. The next step is to evaluate the invariant vector of the 
global propagator with eigenvalue one. Obtaining such a vector is an internal verification 
on the validity of the approximation. The energy E component of the eigenvector at two 
points B and A in Fig. [2] lead to Q c = E B — Ea- 
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1. In the first combination the adiabats will be approximated by Eq. (II 6p . while the hot 
isomagnetic segment will be approximated by Eq. ( 125]) . and the cold isomagnetic by 
Eq. (fTDl . Denoting e~^ TcTc ^ = a and cos(f2 c r c ) = cc, sin(f2 c T c ) = ss, the heat removed 
from the cold bath in this approximation, Q^>pr\ b ecomes: 



Q a PP n ^ ar^ssE^T.-nlE^cc + E^anl) (2ar h 2 r^(q - cc)) 

^ c ~~ (a 2 - 2acc +1) (a - 1) 1 J 

The second term in Eq. (I25p is typically two order of magnitude smaller than the first 
term, therefore we neglect it. The final approximation in this case becomes: 



a PP n b ^ _ arl{n c ssEl"T h - tf h E?cc + jjgajg) 
^ c ~ (a 2 -2acc+l) 



The approximations of Eq. (125]) and (126]) correspond to the data of Fig. [TJ The 
approximation was compared to the numerical calculations in the range of parameters 
corresponding to Fig. HJ in a neighborhood of r c , and T c . A good agreement was 
obtained up to a constant numerical factor. 

The heat removed Q c in Eq. (126|) can change sign meaning that for certain values of 
parameters the refrigeration stops. Fig. shows the alternating cooling as a function 
of t c . The switching points of approximation Eq. (126"]) are conjectured at the points 
where the functions exp(r c r c ) and cos(f2 c T c ) cut each other, (Bottom of Fig. [6]). As 
can be seen these points form a good approximation. 

Fig. HO] is a map of showing regions in parameter space where refrigeration takes place 
and regions where there is no refrigeration. In this map r c = r^, which is a border 
case of the approximation, because the approximation in Eq. ([26]) assumes Tc/t^ 3> 1. 



2. In the second case the adiabats are approximated by Eq. (]T7]) . On both isomagnets a 
second order approximation is used, Eq. (124]) . leading to: 

Qappr2 w (g - ^)(^ - j(r> c ) 2 )(F fe r fe - 0.5(r fe r fe ) 2 ) 

(r c r c + r^) - |(r c r c + r^,) 2 

Fig. |5] shows cycles corresponding to the conditions of Eq. (127]) . with the additional 
condition: u^r^ = constant = 6.252. Comparison of the approximation of Eq. (I27p to 
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FIG. 5: (Color online) Top(a): Cycle trajectories in the H,L,C space which complete approxi- 
mately one revolution on the hot isomagnetic segment, oj^t^ = 6.252. The data for the three 
cycles are: J = 2.5, T h = 10, T c = 9.,uj c = 2.5 k h |= 0.36, k c |= 0.328 r hc = r ch = 0.065625, r c = 
0.008375, iFrom large, to small cycles: uj h = 156.3,78.15,39.075 and therefore r h = 0.04,0.08,0.16 
Bottom(/3): The same cycles in the (Q,Se), (£1,Svn)> planes. 

numerical values of Q c show good agreement with deviations up to ~ 20%. Section IVl 
addresses a large subfamily of cycles corresponding to Eq. (|2T|) . with the additional 
condition of OJhTh = constant, see Fig. [TTJ 

3. In the third approximation we will distinguish between two cases. 
Case A; 

First we combine Eq. ( I18p for the adiabats and Eq. (J22J) for the isomagnets. This case 



becomes equivalent (for (T c r c + T h T h ) < 2) to the frictionless cycles studied before [29 ]. 
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This means that there is no coupling between the energy E and the other variables L 
and C. As a result for short times a bang-bang type solution is optimal 9]. Expanding 
the exponents to second order one gets the following expression: 

(S-^)(r c r c rv ft ) 

~ h (T c T c + T h T h )-(T c r c + T h T h y 1 J 

It it can be shown using Eq. flAJ) that Q^ pr3 > for {T c t c + T h r h ) < 2. 
Case B; 

Combining Eq. ( 120]) for the adiabats and Eq. (122]) for the isomagnets leads to case B. 
Computing the eigenvalues of the corresponding global propagator, the variables Ha 
and Hb separates beautifully from the other operator expectation values. The result: 

(fft - tXrcTjv,) - 2|S(i - A c )A^ 2 (i - C ) 

Qappr3b t? ">> ffc /oq\ 

where c = cos( Jr adi ) and v4 ft / c = exp(— r^/cT^). Fig. plots Eq. ( ]29"]) as a function 
of r a( ij with some internal scaling. The approximation demonstrates the sign changes 
of the heatflow. The cycles of Fig. [7] correspond to Eq. ( 129"j) . In Section |IH Dl we will 
further elaborate on the properties of those cycles and their approximations. 



D. The discontinuous character of the sudden cycle families 

An example of the discontinuous behavior is shown in Fig. [7] corresponding to Eq. (129]) . 
When the time allocation on the adiabats is increased the cycles changed from a concave 
shape to a convex shape. Two of these cycles are shown in Fig. [JJ Additional reduction in 
cycle times leads to Q c < 0, which means that the cycles cease to be refrigerators. Then 
by further reducing the allocated times, the cycles suddenly transformed into a concave 
shape. In addition Q c > 0, refrigerator cycles again. Fig. [7J shows an example for the stated 
behavior of the sudden cycles; there are families of cycles with a small 5 neighborhood, 
beyond which a discontinuity emerges. Changing parameters can lead to another small 
neighborhood. In Fig. [TJthe times on the segments were changed proportionally, decreasing 
the overall cycle time, so that r cyde — » 0. Fig. [S] presents the cycles of Fig. [TJin the H, L, C 
space. 
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FIG. 6: (Color online) Top(a): Eq. (|26p is presented when changing the values of r c , for two 
different values, as denoted on the figure: = 0.00075 red line , = 0.00025 green line. 
Bottom(/3): The functions exp(7 c r c ) and cos(S1 c t c ) are shown as functions of r c . Their crossing 
points compare well to the exact points(Top figure). The other parameters are J = 2,r c h = The = 
0.00035, w c = 0.1, u h = 6, T h = 15, T c = 14, as in Fig. [TJ 

IV. THE COEFFICIENT OF PERFORMANCE (COP) AND THE ENTROPY 
GENERATION (S u ) 

The coefficient of performance (COP) is denned as the heat extracted divided by the 
work input: 

COP w Sr- (30) 



19 



0.8 



0.795 




0.79 



> 



c/3 



0.785 



a 




0.786 



0.784 



> 



0.782 







1 ' 1 1 - 

c „ 


\ 
\ 


1 


^ rT^^I (P)_ 

1 , 1 i c 



8 



10 



12 



FIG. 7: (Color online) Top(a): Se and Svn as a function of Q cycles for three different sets of 
times , r h ,T ch = T hc ,T c = are: 0.000405, 0.4194, 0.029538 green (a), 0.00045, 0.466, 0.03282 or- 
ange (b), 0.00036, 0.3728, 0.026256 non refrigeration cycles, 0.000225, 0.233, 0.01641 ma- 
genta (c) cycles. Bottom(/3): Continuation of Top by lowering the cycle time proportion- 
ally. Th,r c h = Th c ,T c = are: magenta (c) as on Top, 0.0001125, 0.1165, 0.008205 green (d), 
0.00005625, 0.05825, 0.0041025 maroon (e) The parameters for all the cycles are: J = 1.25, T h = 

4,T C = 3.6, u h = 11, a; c = 6.5, 4 = 0.36,4 = 0.0656, see Eq. ([U]) 
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FIG. 8: (Color online)The trajectories in the space H, L, C of the cycles of Fig. [7] where the 
cycles are indicated by letters corresponding to the cycles of Fig. [7J 
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FIG. 9: (Color online) Q c as a function of T a( a with some internal scaling as denoted on the figure 
in order to achieve the exact results by the approximation of Eq. (|29p . which approximate the 
cycles of Fig. [7J with the same parameters. The fineness of the sign change is obvious. 

The entropy generation (S u ) for a cyclic process is generated in the baths: 

« - (| + |) (3D 

An explicit approximations for COP and S u for the case of Eq. ([26]) is now evaluated. This 
requires to compute the work input, W°". Using the notation of Fig. |2j 

W on « {E C -E D ) - Q a c pprlb (32) 
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FIG. 10: (Color online) A map of parameter regions of time allocation leading to refrigeration in blue 
cycles and non refrigeration regions in red squares. The cycle of Fig. 1 is located in the blue island 
in the lower left corner. The general parameters are: J = 2., T c = 14, Th = 15, oj c = 0.1, w/j = 6, 
k| = 0.36, Kc = 0.328. The parameters used fit case 1: Eq. ([26]) . 

Similarly to Eq. the values of Eq and Ejj are evaluated from the limit cycle invariant 
vector X of the cycle propagator at points C and D. 

Denoting again e~^ TcT ^ = a and cos(fi c r c ) = cc, sin(f2 c r c ) = ss, the work done on the 
cycle (W on ) becomes: 

T h E e h q T h {\ - a){a 2 - lace + 1) + r 2 h (E e h q T 2 h a 2 ss 2 (2a - 1) - ^ssT h E^a{\ - a)) 

W° n PS - - --! 

(a - I) (a 2 - lace + 1) + Qlr 2 a 2 (a - cc) - a 2 r/ l r/ t (o; - 2aT h r h + 21^x^00) 

(33) 

when terms with third and larger orders of are neglected. Noticing that the lowest order 
of Th in the expressions for Q c were second order. The invested work becomes: 

W on ~ -r^r ft (34) 

Eq. (134"]) shows that the invested work is on the cold — > hot adiabat and is dissipated almost 
exclusivly on the hot isomagnetic segment. 
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FIG. 11: (Color online) The 1/COP as a function of inverse cooling power, 1/V C . The data 
correspond to the conditions of approximation Qc Ppr2 , Eq. (|27p . The data for both figures are: 
J = 2.5, Th = 10, T c = 9, u c = 2.5, r c = 0.008375, T c h = The = 0.065625 and also the condition 
(cjh Th = 6.252). Besides that = 625.2 is the largest Uh value for both plots. The differences 
between the plots are the k values. For the orange circle figure k\ = 0.328, = 0.36, where three 
points correspond to the three cycles of Fig. whose values are denoted on the figure. The 
largest Th is 0.32. For the red diamond plot = 3.6, and Kc = 3.28. For this figure the largest 
t h is 0.5157905. 

The COP of the cycle is approximated as: 

COP " Th E-r h (^-2acc + l) (35) 



To compute the approximation for S u , according to Eq. (I3"T1) . we use Qc from Eq. 
Qh requires an additional approximation. To first order Qh = — W on , when W on is given 
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FIG. 12: ln("P c ) as a function of J/T c for the three classes of sudden refrigerators. The cycle 
parameters for type 1 correspond to Fig. 1 where T c /T^ was kept constant and T c was varied. The 
cycle parameters of type 2 correspond to Fig. 2 where T c /Th was kept constant and T c was varied. 
The cycle parameters of type 3b correspond to Fig. 7 where T c /Th was kept constant and T c was 
varied. 

by Eq. ( j34|) . therefore S u becomes: 

ar 2 h (n c ssE^T h -nlE^cc + E^anl) , r h E^Y h \ _ t^T„ 



5 (a*-2acc+l)T c + ^)^J ijT" ( } 



Eq. fl36|) shows that most of the entropy production is generated on the hot isomagnetic 
segment. Fig. [TT] present the dependence of the 1/COP on 1/V C , for the subfamily of cycles 
of Eq. (12 7p . The cycles are chosen with the condition u^t„ = constant, for two different 
heat transport coefficients to the bath. This analysis in inspired by the studies of Gordon 



et al. 
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35[ . where a universal behavior of optimal cycles with preassigned cycle times 
was observed. For the sudden cycles a continuous behavior is only local, nevertheless for the 
particular family by choosing uihJh — constant we found a similar behavior. 



V. THE INFLUENCE OF COOLING 



The cooling power V c = Q c /t as a function of the cold bath temperature J[T C , is shown 
in Fig. [T^J A similar plot for the regular cycles can be found in ref. 29]. Common 
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to both cases is a minimum temperature beyond which the cycle cannot cool any longer. 
This minimum temperature is obtained when the occupation of the exited level after the 
demagnetization segment is larger than the equilibrium value at the cold bath. When T c is 
increased the cooling power increases exponentially with the equilibrium occupation at T c . 
This minimum temperature is different for the different classes of sudden refrigerators where 
type 2 outperforms the other classes in both cooling power and minimum temperature. In 
addition a maximum cooling rate at T c > J is observed, for all the three classes. This is in 



contrast to a monotonic decrease in V c in the regular cycles [22j, |29|. The maximum can be 
attributed to the inability to dissipate at the hot bath in a very short time the heat at the 
hot bath. This is a characteristic of the global behavior of the sudden cycles. 

The maximum in the cooling rate is also reflected in the approximation Eq. ( |26|) where 
Q c ~ J. Therefore there exists a positive root of x = (J/2T C ) in the equation: 

d\n{Q a ™ rlb /r) ^ 1 d{Qr rlb ) . ( , 

d(J/2T c ) Q^P rlb d(J/2T c ) [ ' 

The derivation is constrained by the fact that as the temperature of the cold bath is varied, 
also proportionally the temperature of the hot bath changes. Writing T c /Th = C T , therefore 
not only E^ q and T c are dependent on T c , but also E^ q and T h . After quite tedious compu- 
tation, neglecting terms by order of magnitude considerations, and taking into account that 
in the case Q c ~ J one obtains: 

q^SS = 2 ^" 2 ' + 27ii^y " i ^~ 2 ' = ■ (38) 

which leads to: 

j \ max j / j 



2tj iy-2j&.) (39) 

The maximum in Eq. (1391) is quite delicate, this comes about since the constraints on the 
parameters are very restrictive. In addition Eq. f|59|) is independent of the ratio C T , as well 
as other parameters of the hot segment of the cycle. On the other hand, (J/2T c ) max depends 
on all the parameters of the cold segment. This suggests that the the maximum point of V c 
could be found also for the case where the hot bath has constant temperature. 



VI. DISCUSSION 



One can say, that practicing science involves making order in seemingly disorder. The 
ensemble of the sudden cycles is an extreme example of that statement. For example, we 
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saw both in Fig (fT0|) and Fig <^ that changing the time allocations continuously results 
with a large number of discontinuities. Therefore continuity, one of the main aids in re- 
search, doesn't help in the case of sudden cycles. Also, when one changes slightly the bath 
temperatures, or transition probabilities, or the field values might result with P c < 0. One 
can try to optimize P c as a function of time allocation, a standard procedure for regular 
cycles, generally doesn't work for sudden cycles. The result in most cases will be, that 
while continuously changing the time allocations - a jump will occur from refrigerator to 
non refrigerator. 



The refrigerators studied belong to the family of four stroke Otto refrigerators 
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24| . The working fluid composed from an ensemble of spin pairs, which is a simplified model 



of a working medium composed of magnetic salts. For this model the dynamics is described 
by the equation of motion for the thermodynamical observables. The present study focuses 
on refrigerators with cycle times shorter or much shorter than the internal time scale of the 
working fluid. As a result the different segments of the Otto cycle become interconnected. 
This characteristic results in the density operator p deviating from the typical diagonal form 
in the energy representation. As a consequence, the energy entropy was always much larger 
than the Von Neumann entropy. These cycles termed sudden cycles, settle to a limit cycle, 
typically after a large number of iterations. The state of this limit cycle is the eigenvector 
of the cycle's global propagator with eigenvalue one. This property is exploited to study the 
performance of the cycles using a vector space of a closed set of operators which is sufficient 
to represent the density operator of the limit cycle. 

The total cycle is analyzed through segment propagators which map the vector space 
of operators. Our classification scheme of families of sudden cycles is based on analytical 
approximate expressions for the propagators on each segment. These segment propagators 
were then synthesized to global propagators. Having done that, we computed the approxi- 
mate eigenvectors with an eigenvalue one for each limit cycle. These approximate solution 
demonstrate the global property of the cycle. We find a special continuous subset of cycles 



whose behavior is similar to the universal plot of Gordon et al. [34|, |35j. It was shown in the 
study that the entropy production S u is generated on the boundary of the hot bath. Finally, 
the sudden cycles possess a maximum cooling power as a function of J/T c as opposed to the 
regular cycles where we demonstrated a monotonic decrease of V c . 

Finite time thermodynamics has been devoted to the study of systems far from equilib- 
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rium. For example the energy distribution of the working medium was not in the thermal 
Gibbs state. The present study is characterized by coherence, a deviation from a Gibbs 
state because the state does not commute with the stationary Hamiltonian. 
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Appendix A: Thermodynamical Relations. 

In order to fulfill the second law of thermodynamics, the maximal efficiency, r] max , of a 
heat engine with working fluid of two coupled spins(with the corresponding Carnot efficiency 
relation) is Cf. Q: 



For the reverse operation as a refrigerator with the same working fluid, the basic inequality 
must change its direction, and as a consequence from Eq. lAlj. the constraint on the minimum 
cold bath temperature T c will be 



V 



max 





(Al) 




(A2) 



COP 



Carnot 



COP 



(A3) 



T h -T c 
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